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The most complicated phenomena of equilibrium statistics, phase separations and transitions of 
various order and critical phenomena, can clearly and sharply be seen even for small systems in 
the topology of the curvature of the microcanonical entropy Sb(E, N) = ln[W(E, TV)] (Boltzmann's 
principle BP) as function of the conserved energy, particle number etc.. Consequently, BP allows 
to establish the link toward their microscopic origin and the study of the way how interacting many- 
body systems organize into phase-transitions. Also the equilibrium of the largest possible interacting 
many-body systems like self-gravitating systems is described to great extend by the topology of the 
entropy surface Sb{E, N, L) where L is the angular momentum. Conventional (canonical) statistical 
mechanics describes only a small section out of all equilibrium phenomena in nature and only in 
cases where the so-called "thermodynamic limit" applies (homogeneous phases of "infinite" systems 
interacting with short-range interactions). In this paper I present two examples of phase transitions 
of first order which are of fundamental importance: the liquid to gas transition in a small atomic 
cluster and the condensation of a rotating self-gravitating system into single stars or into multi-star 
systems like double stars an rings. Such systems cannot be addressed by ordinary canonical thermo- 
statistics. I also give a geometric illustration how an initially non-equilibrized ensemble approaches 
the microcanonical equilibrium distribution. 

I. INTRODUCTION 

Many theorems we are used to in conventional macroscopic (canonical) thermo-statistics are wrong 
when statistical mechanics addresses small or other non-extensive systems. Here a revision of the funda- 
mentals is demanded. 

Phase separation of normal systems and also in general the equilibrium of closed non-extensive systems 
are not described by the canonical and grand-canonical ensembles. Only the microcanonical ensemble 
based on Boltzmann's entropy Sb = InW, with W(E) = eotr[S(E — H)], the number of classical or 
quantum states, describes correctly the unbiased uniform filling of the energy-shell in phase-space. 

Various ensembles like the (grand-)canonical arc equivalent to the microcanonical ensemble only if the 
system is infinite and homogeneous, i.e. in a pure and homogeneous phase. Only then exists a one to 
one mapping from the conserved mechanical observables as energy E, particle number N, eventually 
angular-momentum L and others to the intensive variables temperature T, chemical potential // and 
eventually rotational frequency u) and so on. Otherwise the (grand-)canonical ensembles do not reflect 
the equilibrized phase-space distribution of a closed ergodic Hamiltonian system ^, ra, 0] see also 
Barre et al [g). (Grand-) canonical potentials are non-analytic at phase transitions whereas SbCE, N,- ■ ■) 
remains multiply differentiable also there. 

This program is far from only academic interest and is deeply demanded in many fields of condensed 
matter. A pseudo Riemannian geometry must replace Ruppeiner's Riemannian geometry of fluctua- 
tions & M. It leads to negative heat capacities at phase-separation which cannot be explained in am 
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canonical formalism and which are well documented experimentally c. 
as also postulated theoretically since long cf. ||, [l7], [l^, 

Besides conventional extensive systems in the thermodynamic limit, this formulation of thermo- 
statistics addresses additionally the following objects: 




it be explained m any 
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• Astro-physical systems with their typical negative heat capacity are clearly outside of any canonical 
approach g §, ^MM- 

• The same is true for the original goal of Thermodynamics, the description of phase separation 0, 
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• and of course small systems like excited nuclei, atomic clusters etc., which are addressed more 
recently, where many experimental results are now cumulating §, [u], £|, ||, [l6| , and which 
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are exotic from the canonical point of view c. f. section Q. 

II. BOLTZMANN'S ENTROPY S B VS. EQUIPARTITION ENTROPY Sequip. 

The microcanonical ensemble is defined by the uniform (completely unbiased) probability distribution 
on the energy manifold £ : 

p = (i) 

tr[6(E-H)] 

tr[A] := J d m q A(q) (2) 
? = x,P (3) 
Boltzmann's entropy is then (Boltzmann' principle): 

SB^M^J^SP-Htf)]} (4) 

That is the definition of entropy we used in all our previous work. 

There may be another entropy: Here we follow the book of Berdichevsky [[HJ for a little while: 

OH Jd^ q 6[E-H(q)}q^ 

< ft a~ > = =T~ — (5) 

% Jd 6N qd[E-H(q)] 

_ jd eN q e[E-H(^)} 

Jd eN qS[E-H(q)] 

= < Qkj— > 7) 
oqk 

• ^equipart. (S) 

where the integral over the 5-function in eq.(o) was transformed into an integral over the (6N — l)-dim. 
energy surface and then for the step from eq.(£j) to (||) the law of Gauss-Ostogradski was used. As eq.@ 
is independent of which d.o.f is taken, this is the equipartition theorem for a finite system. 
One is now tempted to define an equipartition entropy S equ i pa rt. 



dS, 



equipart. ± 
dE ^equipart. 



(9) 



= jd eN qS[E-H(q)} 
Jd eN qe[E-H{q)} 
S eqmpar t. = ln( , n/ * ,v / d 6N qe[E-H(q)]\. (11) 



iV!(27rS) 3Ar 

However, the microcanonical probability P is still (c.f. eqs.(|l| 

S(E - H) 

P = — ir- 12 

trS(E - H) 
Thus : 

Seqmpart. ± -tr[P ln(P)], (13) 

because P = for E ^ H conform to eq.(|]) and Sb- S eqmpar t. has no interpretation within information 
theory: It depends on regions in phase space inaccessible to the system when the conservation of energy 
holds. Therefore, we decided at variance to ref. |5l] to use the original Boltzmann entropy eq.(|J). The 



3 



relative difference between both is oc 1/N. For medium large systems with short range interactions this 
difference is negligible also in comparison to the surface entropy oc iV -1 / 3 which is the main origin of the 
positive curvature discussed in section Thus here the equipartition law (^) holds but the T equipart . 
differs from T = (dSs/dE)- 1 by a term of the order 1/N. 



III. PHASE-TRANSITIONS MICROCANONICALLY AND GRAND-CANONIC ALLY 

For an orientation, how to generalize the concept of phase transitions to microcanonical finite systems 
we start with the Yang-Lee theory for the grand-canonical ensemble of large, homogeneous systems. The 
grand-canonical ensemble for V — > oo: 

Z(T,VL,V) = (r — dN e -[E-»N-TS{E)]/T (w) 

JJo £o 

T/2 /TOO 

= — dedne-^-^-Ts^VT. (15) 
eo JJo 

cons.+lin.+quadr. 



We put the linear term to : If s(e, n) is concave, c.f. figure([l]), then there is a single stationary point 
e s ,n s with 



1 
T 

E. 

T 



dS_ 
~dE 



e s ,n s 

dS 
dN 



and the Laplace transform can be done in Gaussian approximation to get the free energy density: 

-Tln[Z(T,»,V)} 



V 



e s - fm s - Ts s 



Tln(VrAlVr^)) InV 

v + 0{ — ] 



(16) 
(17) 



Consequently, in the fundamental microcanonical ensemble pure phases are indicated by a negative 
maximum curvature, eigenvalue Ai of the Hessian matrix: 



H(e, n) 

eigen-values: Ai > A 2 

Vl 



dr\Be 

d s 

dedn dn 2 



w 

OS 



(18) 



^-eigenvectors : 



Vl,V 2 



direction of order parameter o (19) 



At phase transition, the grand-canonical partition function Z(T,/i) has a zero in the thermodynamic 
limit (Yang-Lee zero Jp2| ) . Here the grand-canonical energy distribution (the kernel of integral (|i"4|)) is 
bimodal c.f. figure (0)7The pure "liquid" phase, at e\ and the pure "gas" phase at coexist at this 
temperature with equal probability. This is of course only possible if s(e,n) has a positive curvature 
between e\ and e^. I.e. the general microcanonical condition for a phase transition of first order is the 
occurrence of a positive maximal curvature Ai. This can now be generalized to finite systems. 



IV. EQUILIBRIUM OF SMALL EXCITED METAL CLUSTERS 



In this section I will shortly discuss the statistical equilibrium of a typical small self-bound system, 
a cluster of a few hundred iVa-atoms interacting by a realistic many-body (metallic) force at constant 
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FIG. 1: Concave entropy with single stationary point o s . In the limit N — > oo this gives a one-to-one mapping 
of E a , N s to T, (j,. The grand-canonical Z(T, fi) has the same analytic properties as the microcanonical W(E, N). 
This is the condition for a pure phase in the conventional Yang- Lee theory j32j. Consequently, we define a pure 
phase microcanonically also for a finite system by a region of negative maximum curvature of S(E, N, • • •). 



pressure of 1 aim. i.e. with the constraint §§t§^ = latm. As function of the excitation energy a typical 
convex intruder shows the signals of a liquid to gas phase-transition of first order, c.f. fig.(||). The 
derivative of the concave hull of S(E) gives the Maxwell construction of the caloric curve T(E). The 
depth As sur j of the intruder (measured from the concave hull) is a measure of the surface tension c.f.Q]. 
Table (Q) compares the transition temperature T tr (inverse slope of the hull) in Kelvin, the latent heat 
per atom (qi a t = 63 — e L ) in e.V., the boiling entropy (the gain of entropy when one atom is evaporated at 
T tr ), the surface entropy (the depth of the intruder As sur t), the average number N sur f of surface atoms 
summed over all clusters, and finally the surface tension a sur f/T tr — s sur fN{j/N aur f per surface atom 
divided by the transition temperature in dimensionless units, with the known bulk values. 

This figure shows clearly the necessary condition for a phase-separation, the bimodality, the appearance 
of two equal minima of the canonical free energy / = e — Ts at ei and ez, which can only appear when 
the entropy s(e) has a positive curvature and correspondingly a negative heat capacity 



/ dS \2 
\dE> 
d 2 S 



(20) 



in between. 





N 


200 


1000 


3000 


bulk 




T tr [K] 


940 


990 


1095 


1156 




qiat [eV] 


0.82 


0.91 


0.94 


0.923 


Na 


Sboil 


10.1 


10.7 


9.9 


9.267 




A$sur f 


0.55 


0.56 


0.45 






-^V surf 


39.94 


98.53 


186.6 


00 




cr/Ttr 


2.75 


5.68 


7.07 


7.41 



TABLE I: Comparison of the parameters of the liquid-gas transition of a cluster of 200, 1000 and 3000 iVa-atoms 
with those of the bulk at an external pressure of latm.. Details c.f. |4J 
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FIG. 2: Entropy per atom s(e) as function of the energy e per atom for a self bound cluster of 1000 sodium atoms. 
The double tangent between ei and ez is the concave hull of s(e). Its derivative gives the Maxwell construction 
in the caloric curve Tie). 



V. EQUILIBRIUM OF SELF-GRAVITATING AND ROTATING SYSTEMS 

Self-gravitation leads to a non-extensive potential energy oc N 2 . No thermodynamic limit exists for S/N 
and no canonical treatment makes sense. At negative total energies these systems have a negative heat 
capacity. This was for a long time considered as an absurd situation within canonical statistical mechanics 
with its thermodynamic "limit" . However, within our geometric theory this is a simple example of the 
pseudo-Riemannian topology of the microcanonical entropy S(E, N) provided non-gravitational physics 
at high densities is excluded. The fundamental importance of self-gravitating systems calls urgently for 
an abolishment of several un-physical axioms used in conventional statistics [E8L p3| . 



To prove the power of our geometric theory of statistical mechanics I show in figure (|3|) the global 
phase diagram of a system of N self-gravitating particles as function of total energy E and given total 
angular-momentum L. As can clearly be seen, our geometric formulation of statistical mechanics gives 
a quite realistic picture of the phase-transition into mono-star or at higher L into multi star configura- 
tions even though this is far beyond any thermodynamic "limit" . In the region of mixed phases one has 
at any E,L point a separation into mono-star, double-star, ring-like, and even more complicated inho- 
mogeneous configurations and the all space-filling gas phase. In this mixed phase the entropy S(E, N) 
has at least one positive eigen-curvature and the system has a negative specific heat (more generally, a 
negative susceptibility), the characteristic signal of a phase transition of first order. More details are in 
refs [^2], ^i], ^5[. Figure (^) gives a few examples of the rich and realistic configurations of the statistical 
equilibrium of such a self-gravitating and rotating many-body system. This demonstrates the superiority 
of our geometric statistics over any conventional canonical theory. These systems are far larger than any 
conventional thermodynamic "limit" can dream of. 
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FIG. 3: Phase diagram of a microcanonical, rotating self-gravitating system [topology of the Hessian, eq.(|T8|)], in 
the (E, L)-plane. The dashed lines E — L = 1 (left) and E = L (right) delimit the region where the Hessian was 
calculated. In the mixed phase the largest curvature is Ai > 




(d) E = -0.42, L = 0.5 



FIG. 4: Color online: Examples of stationary distributions c(x) occurring inside our spherical box. The angular 
momentum is along the vertical axis. Shown are the contour plot and, above it, the density profile: (a) "sin- 
gle star", (b) "double star", (c) "disk", (d) "ring". C.f.: Physical Review Letters-July 15, 2002, cover-page, 
calculations by Votyakov et. al. Jz| 
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VI. GEOMETRIC ILLUSTRATION OF THE APPROACH TO EQUILIBRIUM AND THE 

SECOND LAW 



As Thermodynamics describes a TV-body systems by a few control parameters M only, which are much 
less than the total number 6N of degrees of freedom, it gives only probabilistic answers (the average of 
some observable < O > over the whole ensemble of all systems with the same values of the M control- 
parameters). Therefore, Thermodynamics describes the evolution of the whole ensemble. There is also 
a geometric interpretation of the evolution of a non-equilibrized ensemble to the microcanonical uniform 
filling of the energy-shell in the 6iV-dim. phase space. 

Even though every trajectory spreads over the available phase space and returns after tp i ncarre , 
different points of the manifold have different tp i ncarre which are normally incommensurable. I.e. the 
ensemble spreads irreversibly over the accessible phase-space. 

Due to the redundancy of the information given by the few (M) control-parameters one cannot distin- 
guish the distribution in phase-space from its direct neighborhood. Therefore, in the case of a strongly 
folded (eventually fractal) non-equilibrium phase-space distribution Boltzmann's entropy is the area of 
the closure of the distribution ]34j] . The area of the closure can be calculated by box-counting (34| ^] c.f. 
fig.(^|). Mathematically, the area of the closure is obtained in the limit of box-sides 8 — ► 0. However, by 
several reasons this limit should not be taken in physical applications (see below). This figure illustrates 
also how a non-equilibrium ensemble develops in time, and how for a mixing dynamics it becomes more 
and more dense in the larger available phase-space, so that the area of its closure e s ^> approaches the 
larger area of the new microcanonical ensemble oc (V a + Vb) . This is the geometric meaning of the 
Second Law of Thermodynamics. Due to the inherent redundant information given by Thermodynam- 
ics, represented here by a finite, non-zero, resolution S of the box-counting, this is achieved in a finite 
equilibration time Q . 
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M(t < t ) 



M(t > t ) 



FIG. 5: At time t < to the system is assumed to be equilibrized and be represented by the compact microcanonical 
set in phase-space M(to) = £(V a ), left side. At times to the accessible volume is suddenly enlarged to V a + Vb- The 
system develops into an increasingly folded but non-crossing "spaghetti" -like distribution M(t) in phase-space 
with rising time t > to after opening the volume Vb. The right figure shows only the early form of the distribution. 
At much later times it will become more and more fractal (Gibbs ink-lines) and finally dense in the new larger 
phase space. The grid illustrates the boxes or size 5 6N of the box-counting method. All boxes which overlap with 
A4(t) contribute to the box-counting volume and are shaded gray. Their number is Ng. Then the box-counting 
area is = N$S d where d is the dimension of the microcanonical manifold. 
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VII. CONCLUSION 

The geometric interpretation of classical equilibrium Statistical Mechanics |j6| by Boltzmann's principle 
offers an extension also to the equilibrium of non-extensive systems. In more fundamental, axiomatic 
terms, it opens the application of Thermo-Statistics to "non-simple" systems which are not (homogeneus) 
fluids or in contact with ideal gases. Surprisingly, but also understandably, this is still an open problem 
c.f. ref. J37j page 50 and page 72. 

Because microcanonical Thermodynamics as a macroscopic theory controls the system by a few, usually 
conserved, macroscopic parameters like energy, particle number, etc. without fixing all 6N degrees of free- 
dom, it is an intrinsically probabilistic theory. It describes all systems with the same control-parameters 
simultaneously. If we take this seriously and avoid the so called thermodynamic limit Qiniv^>ao,N/v=p): 
the theory can be applied to small systems but even to the really large, usually inhomogeneous, self- 
gravitating systems, c.f. [^2], ^4|, p5[ . 

Within the new, extended, formalism several principles of traditional Statistical Mechanics turn out 
to be violated and obsolete. E.g. at phase-separation (at negative heat capacity) heat (energy) can 
flow from cold to hot |28| . Or phase-transitions can be classified unambiguously in astonishingly small 
systems. These are by no way exotic and wrong conclusions. On the contrary, many experiments have 
shown their validity. I believe this approach gives a much deeper insight into the way how many-body 
systems organize themselves than any canonical statistics is able to. The thermodynamic limit clouds 
the most interesting region of Thermodynamics, the region of inhomogeneous phase-separation. 

Because of the only one underlying axiom, Boltzmann's principle eq.(|J), the geometric interpreta- 
tion |?6| keeps statistics most close to Mechanics and, therefore, is more transparent. The Second Law 
(AS > 0) can even be shown to be valid in closed, small systems under quite general dynamical condi- 
tions m. 
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